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Vorrichtung zur Bestimmung und/oder Oberwachung des Volumen- 
und/ Oder Massendurchflusses eines Messmediums 



Die Erfindung bezieht sich auf eine Vorrichtung zur Bestimmung und/oder 
Oberwachung des Volumen- und/ Oder Massendurchflusses eines Mess- 
mediums, das eine Rohrleitung in einer Stromungsrichtung durchfliesst, mit 

10 zumindest zwei Ultraschallwandlern, die Ultraschall-Messsignale in die 

Rohrleitung aussenden und aus der Rohrleitung empfangen, und mit einer 
Regel-/Auswerteeinheit, die den Volumen- und/oder den Massendurchfluss 
des Messmediums in der Rohrleitung anhand der Laufeeitdifferenz der 
Ultraschall-Messsignale in Stromungsrichtung und entgegen der Stromungs- 

15 richtung ermittelt. 

Ultraschall-Durchflussmessgerate werden vielfach in der Prozess- und 
Automatisierungstechnik eingesetzt. Sie erlauben es, den Volumen- und/oder 
Massenstrom eines Mediums in einer Rohrleitung beruhrungslos zu 
20 bestimmen. 



Bei dem zuvor genannten Laufzeit-Differenz-Prinzip wird die unterschiedliche 
Laufzeit von Ultraschall-Messsignalen in Stromungsrichtung und entgegen der 
Stromungsrichtung des Mediums ermittelt und ausgewertet. Hierzu werden die 

25 Ultraschall-Messsignale von den Ultraschallwandlern wechselweise in 
Stromungsrichtung und entgegen der Stromungsrichtung des Mediums 
ausgesendet und vom jeweils anderen Ultraschallwandler empfangen. Aus 
der Laufzeitdifferenz der Ultraschall-Messsignale lasst sich die Fliessge- 
schwindigkeit und damit bei bekanntem Durchmesser der Rohrleitung der 

30 Volumendurchfluss bzw. bei bekannter oder gemessener Dichte des Mediums 
der Massendurchfluss bestimmen. 
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Hinsichtlich der Typen von Messgeraten wird unterschieden zwischen 
Ultraschall-Durchflussmessaufnehmern, die in die Rohrleitung eingesetzt 
werden, und Clamp-On Durchflussmessgeraten, bei denen die Ultraschall- 
5 sensoren von aussen an die Rohrleitung mittels eines Spannverschlusses 
angepresst werden. Clamp-On Durchflussmessgerate sind beispielsweise in 
der EP 0 686 255 B1 , der US-PS 4,484,478 oder der US-PS 4,598,593 
beschrieben. 

10 Bei beiden Typen von Ultraschall-Durchflussmessgeraten werden die 

Ultraschall-Messsignale unter einem vorgegebenen Winkel in die Rohrleitung 
bzw. in das Messrohr, in der/ in dem sich das stromende Medium befindet, 
eingestrahlt und/oder empfangen. Urn eine optimale Impedanzanpassung zu 
erreichen, werden die Ultraschall-Messsignale uber einen Vorlaufkorper bzw. 

1 5 einen Koppelkeil in die Rohrleitung eingekoppelt bzw. aus der Rohrleitung 
ausgekoppelt. Hauptbestandteil eines Ultraschallwandlers ist desweiteren 
zumindest ein piezoelektrisches Element, welches die Ultraschall-Messsignale 
erzeugt und/oder empfangt. 

20 Ubiicherweise handelt es sich bei den Ultraschall-Messsignalen, die fur die 
Volumen- bzw. Massendurchflussmessung verwendet werden, urn 
breitbandige Pulse. Es versteht sich von selbst, dass gerade bei kleinen 
Nennweiten der Rohrleitung bez. des Messrohrs der zeitliche Abstand 
zwischen dem ausgesendeten und dem empfangenen Ultraschall-Messsignal 

25 relativ klein ist. Urn hier eine ausreichende Auflosung und damit eine 

zuverlassige Messung durchfuhren zu konnen, wird das Messsignal mit einer 
Abtastrate abgetastet, die einerseits kleiner ist als die Zeitdauer die zwischen 
dem Aussenden und Empfangen eines Ultraschall-Messsignals liegt und die 
andererseits so klein ist, dass innerhalb der Messimpulsdauer mehrere 

30 Abtastwerte abgetastet werden. Die Abtast- bzw. Samplerate ist folglich relativ 
hoch. Die Abtastwerte bzw. die abgetasteten Amplitudenwerte des Ultraschall- 
Messsignals werden einem A-/D-Wandler zugefuhrt. Eine Regel- 
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/Auswerteeinheit, z.B. eine DSP, verwendet die Samplewerte bzw. die 
Abtastwerte dazu, das empfangene Messsignal durch eine kontinuierliche 
Funktion zu interpolieren bzw. moglichst wirklichkeitsgetreu zu rekonstruieren. 
Mathematisch lasst sich dieser Sachverhalt durch die stetige Funktion f(t) = 
5 f(n«T) = a„ darstellen, wobei n = 1, 2, 3, ... - also eine natQriiche Zahl - ist und 
wobei die Koeffizienten^die zu den Zeitpunkten (n»T) gemessenen 
Amplitudenwerte des Ultraschall-Messsignals darstellen. 

Im einfachsten Fall handelt es sich bei der Funktion um die sukzessive, 
10 lineare Verbindung von jeweils zwei aufeinanderfolgenden Samplewerteh / 
Abtastwerten. Da dieses Verfahren fur Messungen mit einer gehobenen 
Messgenauigkeit im Bereich der Ultraschall-Durchflussmessung nicht 
ausreichend ist, ist es bekannt geworden, die Lagrange-lnterpolation oder die 
noch komplexere Interpolation nach Levenberg-Markart zur Rekonstruktion 
15 des empfangenen Messsignals heranzuziehen. 

Die beste und wohl exakteste theoretische Interpolationsmethode fur aus 
Abtastwerten rekonstruierte Signale besteht fraglos in der Verwendung des 
Shannon-Nyquist Theorems, wonach sich eine begrenzte kontinuierliche 
20 Funktion, z.B. ein Ultraschall-Messimpuls, durch eine unendliche Summe von 
gewichteten Spaltfunktionen sin (x) / x darstellen lasst. Die korrekte Formel 
lautet wie folgt, wobei t als Abkurzung fur die Zeit steht und T die Zeitdauer 
zwischen zwei Abtastwerten reprasentiert: 



Problematisch bei der Anwendung dieser Formel ist, dass die Rekonstruktion 
des Messsignals nur dann hundertprozentig korrekt ist, wenn die Anzahl der 
Abtastwerte unendlich ist. Um das Shannon-Nyquist Theorem in der Praxis 



25 




*(--*) 
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anwenden zu k6nnen, 1st es naturlich erforderlich, die Anzahl der gesampelten 
Messwerte nach oben zu beschranken. Hierbei ist darauf zu achten, dass die 
Anzahl der gesampelten Messwerte gross genug ist, urn eine ausreichende 
und hinreichende Rekonstruktion des empfangenen Messsignals erreichen zu 
konnen. Problem hierbei ist eine moglichst kurze Rechenzeit bzw. eine 
hinlanglich grosse Kapazitat des hinter der Regel-/Auswerteeinheit stehenden 
Mikroprozessors. Bei den heute bekannten Auswerteverfahren fur hoch- 
genaue Durchflussmessungen mittels Ultraschall ist der Energiebedarf fur die 
Lieferung von ('Nahezu'-) Echtzeitmessungen so gross, dass hier nur 
Vierdraht-Messgerate geeignet scheinen. Der Einsatzvon sog. Niedrig- 
Energiegeraten, insbesondere von Zweileiter-Ultraschall-Durchflussmess- 
geraten, war bislang aufgrund des hohen Energiebedarfs nicht moglich. Wie 
bereits angedeutet, ist der hohe Energiebedarf vorrangig durch die grosse 
Rechnerkapazitat des Mikroprozessors bzw. der DSP gefordert. Letztlich ist 
der hohe Energiebedarf eine Folge der aufwendigen Auswerteverfahren, 
welche fQr hochdynamische Messungen - insbesondere im Bereich der 
Echtzeit-Messung - mit hoher Messgenauigkeit erforderlich sind. 

Der Erfindung liegt die Aufgabe zugrunde, ein Ultraschall-Durchflussmess- 
gerat mit niedrigem Energieverbrauch vorzuschlagen, 

Die Aufgabe wird dadurch gelost, dass die Regel-/Auswerteeinheit mehrere 
Abtastwerte («,mit i = 1, 2, 3, ...) eines empfangenen Messsignals in 
definierten Zeitpunkten (t) eines vorgegebenen Zeitbereichs ermittelt, 
dass die Regel-/Auswerteeinheit den vorgegebenen Zeitbereich des 
Messsignals durch eine stetige Funktion (f(t)) interpoliert, wobei die stetige 
Funktion (f(t)) durch eine Summe einer vorgegebenen Anzahl (n e N) von 
Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines 

Abtastwerts mit einer Spaltfunktion (^1) mit einer Gauss'schen 

X 

Glockenkurve (e-"* 1 , a e R) entspricht. Wie bereits an vorhergehender Stelle 
erlautert, besteht das Problem mit der Funktion sin (x) / x darin, dass sie fur 
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praktische Anwendungen viel zu langsam gegen Null konvergiert. Die 
Funktion sin(x)/x ist quasi erst bei Minus-Unendlich bzw. bei Plus-Unendlich 
gleich Null. Durch die Produktbildung mit der schnell abfallenden 
zuvorgenannten Gauss'schen Glockenkurve kann dieses Problem elegant 
umschifft werden, wodurch die Rechenkapazitat der Regel-/Auswerteeinheit 
ebenso wie deren Energiebedarf erheblich reduziert werden kann. Die 
Einsparung liegt im Hinblick auf die bislang bekannt gewordenen Verfahren 
bei einem Faktor, der sich in der Grossenordnung von 10 bis 100 bewegt. 

Gemass einer vorteilhaften Weiterbildung der erfindungsgemassen 
Vorrichtung ist vorgesehen, dass die Regel-/Auswerteeinheit zwischen den 
Abtastwerten zumindest einen zusatzlichen Abtastwert bestimmt und diesen 
Abtastwerte bzw. diese Abtastwerte durch die stetige Funktion approximiert, 
wobei die stetige Funktion durch die Summe einer vorgegebenen Anzahl (n e 
N) von Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt 

eines Abtastwerts mit einer Spaltfunktion {^^-) mit einer Gauss'schen 

Glockenkurve (e"^ ,a e R) entspricht. Dieses Verfahren ist in anderem 
Zusammenhang bereits bekannt und tragt den Namen: Oversampling. Durch 
die Zwischenabtastung und Zwischenwertberechnung lasst sich eine bessere 
Auflosung des abgetasteten empfangenen Messsignals und damit eine 
hohere Messgenauigkeit bei der Bestimmung des Volumen- oder Massen- 
durchflusses erreichen. 

Weiterhin schlagt eine gunstige Ausgestaltung der erfindungsgemassen 
Vorrichtung vor, dass die Regel-/Auswerteeinheit einen Abszissenwert (t) 
bestimmt, bei dem ein Ordinatenwert der stetigen Funktion (f(t)) einen 
vorgegebenen Grenzwert erreicht. Bevorzugt handelt es sich bei dem 
vorgegebenen Grenzwert der stetigen Funktion (f(t)) urn ein Maximum. Es 
kann sich jedoch ebenso urn einen Nullpunkt, ein Minimum oder urn einen 
Wendepunkt handeln. 
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Bevorzugt ermittelt die Regel-/Auswerteeinheit anhand der ersten Ableitung 
f (t) der stetigen Funktion f(t) den Abszissenwert (tmax, tmin) im Maximum 
und/oder im Minimum. Dieser ermittelte Abszissenwert ist nachfolgend jeweils 
der direkte Bezugswert fur die Laufzeit des Messsignals in Stromungsrichtung 
bzw. entgegen der StrSmungsrichtung. 

Eine besonders vorteilhafte Ausgestaltung der erfindungsgemasseh 
Vorrichtung sieht vor, dass die Regel-/Auswerteeinheit den Abszissenwert 
(tmax), bei dem die stetige Funktion ein Maximum erreicht, durch eine lineare 
Interpolation der ersten Ableitung der stetigen Funktion (f(t)) nach folgender 
Formel erhalt, und wobei to den Abszissenwert einer ersten Schatzung 
kennzeichnet, bei dem im Zeitintervall (to - T, to + T) ein Maximum oder 
Minimum gemessen wird, und wobei f"(t) die zweite Ableitung der stetigen 
Funktion (f(t)) reprasentiert. Mathematisch lasst sich dieser Sachverhalt durch 
die folgende Formel ausdrucken: 



f\tO) 



Im Zusammenhang mit der erfindungsgemassen Vorrichtung wird es als 
besonders giinstig angesehen, wenn die Regel-/Auswerteeinheit zwei 
Ultrasschall-Messsignale in zwei Zeitbereichen miteinander korreliert, die 
entsprechende diskrete Abtastung von Korrelationspunkten durch eine stetige 
Funktion (f(t)) interpoliert und den Abszissenwert der stetigen Funktion (f(t)) 
bestimmt, bei dem der Ordinatenwert das Maximum erreicht, wobei der 
Abszissenwert ein Mass fur die Zeitverschiebung zwischen den in 
Stromungsrichtung und entgegen der StrSmungsrichtung gesendeten und 
empfangenen Ultraschall-Messsignale ist. Das Ergebnis dieser Auswertung 
liefert somit direkt die Zeitdifferenz zwischen den beiden in unterschiedliche 
Richtungen ausgesendeten und empfangenen Ultraschall-Messsignalen. 
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Die Messgenauigkeit des Ultraschall-Durchflussmessgerats hangt ent- 
scheidend von der korrekten bzw. optimalen Wahl des Koeffizienten a der 
Gauss'schen Glockenkurve ab. Urn zu gewahrleisten, dass der Koeffizient 
a optimal bestimmt ist, wird er gemass einer vorteilhaften Ausgestaltung der 
erfindungsgemassen Vorrichtung in Abhangigkeit von der Anzahl der 
Messpunkte (MaxSamplei) bestimmt. Insbesondere ermittelt die Rechen- 
/Regeleinheit jeweils einen optimalen Wert fur den Koeffizienten a in 
Abhangigkeit von der Anzahl der Messpunkte (MaxSamplei) mittels eines 
mathematischen Simulationsprogramms. 

Bevorzugt ist eine Speichereinheit vorgesehen ist, in der jeweils der optimale 
Wert fur den Koeffizienten a in Abhangigkeit von der Anzahl der Messpunkte 
(MaxSample) abgespeichert ist. Beispielsweise sind die Werte in einer Tabelle 
abgelegt. Mit dieser Methode lasst sich wiederum Rechenzeit und Energie 
sparen, da der jeweilige Wert bei Bedarf einfach aus der Tabelle entnommen 
werden kann. In Abhangigkeit von der jeweils geforderten Messgenauigkeit 
und/oder in Abhangigkeit von der-momentan zur Verfugung stehenden 
Energie lasst sich so stets ein optimales Messergebnis erzielen. 

Aufgrund des geringen Energieverbrauchs ist es moglich, das erfindungs- 
gemasse Ultraschall-Durchflussmessgerat als Zweileiter-Durchflussmessgerat 
auszubilden. Zweileiter-Technologie bedeutet, dass die Energieversorgung 
des Gerats und die MesswertQbermittlung an eine entfernte Leitstelle und ggf. 
die Konfigurierung und Parametrisierung des Gerats von der entfernten 
Leitstelle her uber lediglich zwei Leitungen erfolgt. Da die Verdrahtungskosten 
Qblicherweise einen relativ hohen Anteil der Gesamtkosten ausmachen, 
lassen sich hier beachtliche Einsparungen erzielen. Aufgrund des geringen 
Energieverbrauchs ist es dariiber hinaus auch moglich, das 
Durchflussmessgerat mit einer internen Energiequelle auszustatten. Die 
Kommunikation mit einer entfernten Leitstelle kann dann uber Leitungen oder 
drahtlos, z.B. Qber Funk, erfolgen. 
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Die Erfindung wird anhand der nachfolgenden Figuren naher erlaiitert: Es 
zeigt: 

Fig. 1: eine schematische Darstellung eines Diagramms, bei dem die 
abgetasteten Amplitudenwerte des Messsignals gegen die Zeit aufgetragen 
sind (-> Stand der Technik), 

Fig. 2: eine schematische Darstellung eines Diagramms, an dem die 
Arbeitsweise der erfindungsgemassen Vorrichtung visualisiert wird und bei 
dem die abgetasteten Amplitudenwerte des Messsignals gegen die Zeit 
aufgetragen sind und 

Fig. 3: eine schematische Darstellung einer bevorzugten Ausgestaltung der 
erfindungsgemassen Vorrichtung. 

In Fig. 1 sind die Amplitudenwerte a„ von sechs Abtastwerten, die im 
Zeitabstand T abgetastet wurden, gegen die Zeit aufgetragen. Ober die 
gestrichelte Linie ist die lineare Verbindung zwischen zwei Messpunkten 
dargestellt. Die durchgezogene Kurve f(t) kennzeichnet eine Kurve, die nach 
einem der bekannten Rekonstruktions-Algorithmen, z.B. nach Lagrange oder 
nach Levenberg-Markart gefunden wurde. Nachteil der bekannt gewordenen 
Methoden ist der hohe Rechenaufwand und damit einhergehend der hohe 
Energieverbrauch der Regel-VAuswerteeinheit bzw. des Mikroprozessors oder 
der DSP. 

Erfindungsgemass ermittelt die Regel-/Auswerteeinheit 11 mehrere Abtast- 
werte (o,mit i = 1, 2, 3, ...) eines empfangenen Messsignals in definierten 
Zeitpunkten (t) eines vorgegebenen Zeitbereichs. Anschliessend werden die 
Abtastwerte in dem vorgegebenen Zeitbereich durch eine stetige Funktion 
(f(t)) interpoliert, wobei die stetige Funktion (f(t)) durch eine Summe einer 
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vorgegebenen Anzahl (n e N) von Wavelets (W) gebildet ist und wobei jedes 
Wavelet (W) dem Produkt eines Abtastwerts mit einer Spaltfunktion (^^1) 

mit einer Gauss'schen Glockenkurve {e°* , a e R) entspricht. 
Im Zeitbereich hat die Formel die folgende Form: 

- r r 1 ' « sin[^-H)l r^f 
»=— L 1 J »=-~ 7t(—-ri) 

Mittels der zuvorgenannten Berechnung lasst sich eine effiziente und schnelle 
Interpolation bei relativ geringer Rechnerkapazitat erzielen. 

Betrachten wir den Fall des Oversampling und nehmen an, das die 
Oversampling-Rate gleich a ist, wobei a eine ganze Zahl und grdsser Oder 
gleich 2 ist. Folglich ist jeder Abtastbereich der Zeitdauer T in a Teil- 
Abtastbereiche der Zeitdauer T/o unterteilt. 

Im folgenden konzentrieren wir uns auf den Zeitbereich -TV2 < t < +T/2. 
Selbstverstandlich konnen auch andere Werte von t analog berechnet werden 
indem der Abtastbereich entsprechend verschoben wird. Die Formel fur die 
Funktion f(iT/o) lautet fur i im Bereich: -(a/2+1) < i <> +(a/2), wobei i * 0 sein 
muss: 



sin 

/£>-£< 



a £L K( L_ n) — a 

a 

c(i/o\n) konvergiert mit wachsendem n sehr schnell gegen Null, so dass in der 
Praxis die Naherung gemacht werden kann, dass c(i/a, n) gleich Null ist. 
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sobald der absolute Wert von n einen vorgegebenen Grenzwert erreicht. 
Dieser Grenzwert hangt letztlich von der geforderten bzw. notwendigen 
Messgenauigkeit ab und wird nachfolgend als MaxSamples gekennzeichnet. 
Es hat sich gezeigt, dass eine gewlinschte Messgenauigkeit sich in der Praxis 
erreichen lasst, wenn der Wert von MaxSamples im Bereich von ca. 3 bis 10 
liegt. 

Damit ergibt sich die folgende Gleichung: 

•rp ^-MaxSamples 

/(§)= X<v c (;;>"> 

Die Koeffizienten c(i/o,n) werden einmal berechnet und dann in einer Tabelle 
der Gr6sse (2* MaxSamples + 1) * ( o - 1) abgespeichert. Somit lasst sich 
nahezu einer Verarbeitung der Messwerte in Echtzeit erreichen. 



n 


i/c = -0.25 


i/c = 0.25 


i/o = 0.5 


-5 


XI 


X2 


X3 


















0 


Yl 


Y2 


Y3 


















5 


Zl 


Z2 


Z3 



Wesentlich ist es, anhand der in einem Zeitbereich T abgetasteten 
Amplitudenwerte den Zeitpunkt t max zu ermitteln, an dem ein relativer 
maximaler Amplitudenwert (oder auch ein relativer minimaler Amplitudenwert) 
der Funktion f(t) auftritt. Beispielsweise kann dieser Wert zwischen zwei 
aufeinanderfolgenden Abtastwerten liegen. Ist Wbekannt, so ist es auch 
moglich das tatsachliche Maximum f(w) (oder Minimum) der Funktion f(t) zu 
bestimmen. 
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Die Vorgehensweise ist die folgende (siehe auch Fig. 2): Zuerst wird ein 
Maximum (oder Minimum) in dem Zeitbereich bestimmt, indem die 
Abtastwerte miteinander verglichen werden. Nehmen wir an, dass ao dem 
maximalen Wert (oder dem minimalen Wert) entspricht. Zudem wird 
5 angenommen, dass mindestens MaxSamples vor und hinter ao vorhanden 
sind. Fur praktische Anwendungen liegt der Parameter MaxSamples im 
Bereich von 3 -10. 



10 Bei dem Abszissenwert t m ax ist die erste Ableitung f (t) der Funktion f(t) also 
gleich Null. Die Formel lautet demnach: f(W) = 0. Somit ist a 0 das Maximum 
(oder das Minimum) von dem Zeitbereich bzw. Zeitintervall [ad, wobei gilt: -T 
< t max <T. 

Der Wert t max ISsst sich durch eine lineare Interpolation der ersten Ableitung 
1 5 gemass der folgenden Formel finden: 

max y,,^ 

Wie bereits gesagt wird erfindungsgemass folgendes Interpolations-Wavelet 
20 genutzt: 

. . sin(x) _„> 
g(x)-— —e 
x 

Die erste Ableitung g'(x) und die zweite Ableitung g"(x) der Funktion g(x) 
25 lauten: 

ff M W = r2^-2^ + (2«-l)^-4«.cos(x) + 4^.sin(x)V 

30 

Eine Grenzwertberechnung mit x = 0 fuhrt zu g(0) = 1 , g'(0) = 0, 
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10 



g"(0) = (-2cc-1)/3 

Hieraus ergeben sich fOr f(t), f(t), f'(t) folgende mathematischen Gleichungen: 
Folglich gilt fur t = 0 



15 



/"(0)=- F 



/*(0)=-^I(-i) n -^^ v 

a 0 (2a+ \)7t 2 +2X(-l)"« n -(2«^ 2 +^)e 



, = A"). 



«=1 " 



"max /■!! 



'"<°> o.(2a+V +22(-l)"(«_„ + a n ).(2fl*r 2 +-j)e^ 

•3 «=1 



20 Wird weiterhin angenommen, dass nur die Abtastwerte innerhalb des Sub- 
Zeitbereichs von a. Ma xsampies bis a Ma xsampi e s relevant sind, dann wird a n =0, 
wenn |n| > MaxSamples ist. Die letzte Gleichung lautet dann: 



25 



'maxl "" T 



Ma xSam ples 



a- n -a„ c wprV 



i MaxSamples 1 2 2 

a 0 (2«+> 2 +2 X (-D-Ca^+aJ^^-K^" 



Wenn der Bereich [ad in den Bereich [a r a 0 ] vertikal verschoben wird, so lasst 
sich die zuletzt genannte Formel vereinfachen und man erhait: 
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max 2 MaxSamples 1 2 2 

2 £ (-D"(^+^-2«oM2«^ 2 +— 



Um den optimalen Wert fur den Koeffizienten a zu finden, sollte Wi = t m ax2 
gesetzt werden, da sich der Abszissenwert des Maximums bei einer 
translatorischen Verschiebung nicht andern sollte, da er gegeniiber einer 
Translation invariant sein musste. Allerdings ist in der strengen Theorie t max i * 
tmax2 (Qbrigens sogar dann, wenn MaxSamples = - ist). Grund ist, dass die 
Funktion, die den verschobenen Bereich [a r a 0 ] interpoliert keine einfache 
Translation der Funktion ist, die den Bereich [aj interpoliert. Daher ist die 
Annahme, dass t max i * W ist, nur giiltig fur das nachfolgend formulierte 
Residuum: 

| MaxSamples | ^ 2 2 

residue(a,MaxSamp!es) = 2a+- + 4 J) (-l) n (2a+- I - r K <w " 

Der absolute Wert von diesem Residuum ist ein Zeichen dafUr, wie gut das 
gewahlte Wavelet fur den praktischen Gebrauch ist. Die Analyse des 
Residuums mittels eines Mathematischen Simulationsprogramms (z.B. 
Mathcad) ermoglicht es, einen optimalen Wert fur a fur eine vorgegebene 
Anzahl von (2xMaxSamples) um aoherum zu finden. Optimale Werte fur liegen 
in der GrSssenordnung von 0.01 bis ca. 0.04. Wird ein Wert fOr a eingesetzt, 
der abseits des optimalen Wertes liegt, so ist die Interpolationsfunktion 
zwischen zwei Abtastwerten nicht 'glatt', bzw. sie hat signifikante Harmonische 
jenseits der Nyquist Grenze, aber geringe Abweichungen vom optimalen Wert 
(z. B. 5%) haben wenig Einfluss auf das Ergebnis. Letzteres zeigt wie gut die 
vorgeschlagene Vorrichtung zur Durchfuhrung der praktischen Messungen 
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geeignet ist. Insbesondere beeinflussen Rundungsfehler das Messergebnis 
minimal. 

In Praxis werden zwei Koefflzienten-Tabellen [c{\ and [di] (0 < i< MaxSamples), 
aufgebaut. Dann ergibt sich folgender Zusammenhang: 



MaxSamples 

t = T ~ mit 

max * MaxSamples 

c ™d </„ =-2.(-l)".(2ow 2 + A.) e - flwl( " J - I) 

Die Koeffizientenbereich [q] and [di] weden normiert, so dass gilt: ci=1 . Der 
optimale Wert fur a, der von der Anzahl der Abtastwerte MaxSamples 
abhangt, wird in eine Tabelle eingelesen. Die Koeffizienten werden berechnet 
und in einer Tabelle abgespeichert. Im Betrieb bedeutet das fur den Mikro- 
prozessor bzw. die Regel-/Auswerteeinheit des Ultraschall-Durchflussmess- 
gerats, dass nur einfache Rechenoperationen wie Additionen, 
2xMaxSamples Multiplikationen und eine Division durchzufuhren sind. 

Wie bereits an vorhergehender Stelle mehrfach erwahnt, wird zur Messung 
des Durchflusses mittels eines Ultraschall-Durchflussmessgerats ein Laufzeit- 
differenzverfahren angewendet. Ein Ultraschall-lmpuls wird in Stromungs- 
richtung (Up) des Messmediums in die Rohrleitung bzw. in das Messrohr 
eingestrahlt, von einem Ultraschallwandler empfangen und nachfolgend 
werden von einem (high-speed) schnellen A-/D-Wandler innerhalb eines 
vorgegebenen Zeitbereichs mehrere Abtastwerte [upj gesammelt. Dasselbe 
Signal wird dann entgegen der Stromungsrichtung (Down) in die Rohrleitung 
bzw. in das Messrohr ausgesendet, ebenfalls von einem Ultraschallwandler 
empfangen und von dem A-/D-Wandler gesampelt. Innerhalb eines 
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vorgegebenen Zeitbereichs werden ebenfalls mehrere Abtastwerte [dnj] 
gesammelt. Die Zeitdifferenz zwischen den beiden Messsignalen ist 
proportional zur Fliessgeschwindigkeit des Messmediums in der Rohrleitung. 
Die beiden Abtastwerte werden miteinander korreliert gemass der Formel: 

5 

cotTi = Yup„dn^, 

Der maximale Wert dieses Bereichs korrespondiert zu der Zeitdifferenz der 
10 beiden zuvorgenannten Ultraschall-Messsignale. Fur eine hochgenaue 

Berechnung des Volumenstsroms ist dieses Korrelationsverfahren jedoch viel 
zu ungenau. Daher wird auf die Korrelationsfunktion das zuvorgenannte 
erfindungsgemasse Verfahren angewendet. 

15 Fig. 3 zeigt eine schematische Darstellung der erfindungsgemassen 

Vorrichtung in der Ausgestaitung als Inline-Ultraschall-Durchfiussmessgerat 1. 
Das Ultraschall-Durchflussmessgerat 1 ermittelt den Volumendurchfluss bzw. 
den Massendurchfluss des Messmediums 4, das in Stromungsrichtung (S 
bzw. Up) in der Rohrleitung 2 fliesst, nach der bekannten Laufzeitdifferenz- 

20 Methode. 

Wesentliche Komponenten des Inline Ultraschall-Durchflussmesserats 1 sind 
die beiden Ultraschallwandler 5, 6 und die Regel-/Auswerteeinheit 11. Die 
beiden Ultraschallsensoren 5, 6 sind mittels einer in der Fig. 1 nicht gesondert 
25 dargestellten Befestigungsvorrichtung in einem Abstand L voneinander an der 
Rohrleitung 2 angebracht. Entsprechende Befestigungsvorrichtungen sind aus 
dem Stand der Technik hinlanglich bekannt und werden auch von der 
Anmelderin angeboten und vertrieben. Die Rohrleitung 2 hat einen vorge- 
gebenen Innendurchmesser di . 

30 

Ein Ultraschallwandler 5; 6 weist als wesentlichen Bestandteil zumindest ein 
piezoelektrisches Element 9; 10 auf, das die Ultraschall-Messsignale erzeugt 



FL0226-DE 



und/oder empfangt. Die Ultraschall-Messsignale werden jeweils uber die 
Koppelelemente 7, 8 der beiden Ultraschallwandler 5; 6 in die vom Medium 4 
durchstromte Rohrieitung 2 eingekoppelt bzw. aus der vom Medium 
durchstromten Rohrieitung 2 ausgekoppelt. Das Koppelelement 7, 8 sorgt fur 
eine mSglichst gute Impedanzanpassung der Ultraschall-Messsignale beim 
Ubergang von einem Medium in das andere. Mit SP ist Qbrigens der 
Schallpfad gekennzeichnet, auf dem sich die Ultraschall-Messsignale in der 
Rohrieitung 2 bzw. in dem Messmedium 4 ausbreiten. Im gezeigten Fall 
handelt es sich urn eine sog. Ein-Traversenanordnung, in der die Ultraschall- 
wandler 5, 6 angeordnet sind. Eine Traverse kennzeichnet dabei den 
Teilbereich des Schallpfades SP, auf dem ein Ultraschall-Messsignal die 
Rohrieitung 2 einmal quert. Die Traversen konnen je nach Anordnung der 
Ultraschallwandler 5, 6 und ggf. unter Einfiigen eines Reflektorelements in 
den Schallpfad SP diametral oder chordial in der Rohrieitung bzw. in dem 
Messrohr 2 verlaufen. 
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Patentanspruche 

5 

1 . Vorrichtung zur Bestimmung und/oder Oberwachung des Volumen- und/ 
Oder Massendurchflusses eines Messmediums (4), das eine Rohrleitung (2) in 
einer Stromungsrichtung (s; Up) durchfliesst, mit zumindest zwei Ultraschall- 
wandlern (5, 6), die Ultraschail-Messsignale in die Rohrleitung (2) aussenden 

10 und empfangen, und mit einer Regel-/Auswerteeinheit (11), die den Volumen- 
und/oder den Massendurchfluss des Messmediums (4) in der Rohrleitung (2) 
anhand der Laufeeitdifferenz der Ultraschail-Messsignale in Stromungs- 
richtung (S; Up) und entgegen der Stromungsrichtung (Down) ermittelt, 
dadurch gekennzeichnet, 

15 dass die Regel-/Auswerteeinheit (11) mehrere Abtastwerte (a,. mit i = 1, 2, 3, . 
...) eines empfangenen Messsignals in definierten Zeitpunkten (t) eines 
vorgegebenen Zeitbereichs ermittelt, 

dass die Regel-/Auswerteeinheit (11) den vorgegebenen Zeitbereich des 
Messsignals durch eine stetige Funktion (f(t)) interpoliert, wobei die stetige 
20 Funktion (f(t)) durch eine Summe einer vorgegebenen Anzahl (n e N) von 
Wavelets (W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines 

Abtastwerts mit einer Spaltfunktion ( sin ^ ) mit einer Gauss'schen 

x 

Glockenkurve (e" 0 * 2 , a e R) entspricht. 

25 2. Vorrichtung nach Anspruch 1 , 
dadurch gekennzeichnet, 

dass die Regel-/Auswerteeinheit (11) zwischen den Abtastwerten zumindest 
einen zusatzlichen Abtastwert bestimmt und diesen Abtastwerte bzw. diese 
Abtastwerte durch die stetige Funktion approximiert, wobei die stetige 
30 Funktion durch die Summe einer vorgegebenen Anzahl (n e N) von Wavelets 
(W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Abtastwerts 
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mit einer Spaltfunktion ( sm ^ ) mit einer Gauss'schen Glockenkurve (e"** 2 ,a 

x 

e R)entspricht 

3. Vorrichtung nach Anspruch 1 Oder 2, 
5 dadurch gekennzeichnet, 

dass die Regel-/Auswerteeinheit (11) einen Abszissenwert (t) bestimmt, bei 
dem ein Ordinatenwert der stetigen Funktion (f(t)) einen vorgegebenen 
Grenzwert erreicht. 

10 4. Vorrichtung nach Anspruch 3 
dadurch gekennzeichnet, 

dass es sich bei dem vorgegebenen Grenzwert der stetigen Funktion (f(t)) urn 
einen Nullpunkt, ein Maximum, ein Minimum oder urn einen Wendepunkt 
handelt. 

15 

5. Vorrichtung nach Anspruch 1, 
dadurch gekennzeichnet, 

dass die Regel-/Auswerteeinheit (11) anhand der ersten Ableitung f (t) der 
stetigen Funktion f(t) den Abszissenwert (tmax, tmin) im Maximum und/oder 
20 im Minimum bestimmt. 

6. Vorrichtung nach Anspruch 1, 
dadurch gekennzeichnet, 

dass die Regel-/Auswerteeinheit (11) den Abszissenwert (tmax), bei dem die 
25 stetige Funktion ein Maximum erreicht, durch eine lineare Interpolation der 

ersten Ableitung der stetigen Funktion (f(t)) nach folgender Formel erhalt, und 
wobei tO den Abszissenwert einer ersten Schatzung kennzeichnet, bei dem im 
Zeitintervall (tO - T, tO + T) ein Maximum oder Minimum gemessen wird, und 
wobei f '(t) die zweite Ableitung der stetigen Funktion (f(t)) reprasentiert: 
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f'itO) 



7. Vorrichtung nach einem oder mehreren der vorhergehenden Anspruche, 
dadurch gekennzeichnet, 

dass die Regel-/Auswerteeinheit (11) zwei Ultraschall-Messsignale in zwei 
Zeitbereichen miteinander korreliert, die entsprechende diskrete Sammlung 
von Korrelationspunkten durch eine stetige Funktion (f(t)) interpoliert und den 
Abszissenwert der stetigen Funktion (f(t)) bestimmt, bei dem der 
Ordinatenwert den Maximalwert erreicht, wobei der Abszissenwert ein Mass 
fur die Zeitverschiebung zwischen den in Strfimungsrichtung (S, Up) und 
entgegen der Stromungsrichtung (Down) gesendeten und empfangenen 
Ultraschall-Messsignale ist. 

8. Vorrichtung nach Anspruch 1 oder 2, 
dadurch gekennzeichnet, 

dass die Rechen-/Regeleinheit (1 1) jeweils einen optimalen Wert fUr den 
Koeffizienten (a ) in Abhangigkeit von der Anzahl der Messpunkte 
(MaxSamplei) mittels eines mathematischen Simulationsprogramms 
bestimmt. 

9. Vorrichtung nach Anspruch 8, 
dadurch gekennzeichnet, 

eine Speichereinheit (12) vorgesehen ist, in der jeweils der optimale Wert fur 
den Koeffizienten (a ) in Abhangigkeit von der Anzahl der Messpunkte 
(MaxSample) abgespeichert ist. 
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Zusammenfassung 

Die Erfindung beziehtsich auf ein Ultraschall-Durchflussmessgerat (1), das 
sich durch einen niedrigen Energieverbrauch auszeichnet. 

Erfindungsgemass ermittelt die Regel-/Auswerteeinheit (11) mehrere 
Abtastwerte (a,mit i = 1, 2, 3, ...) eines empfangenen Messsignals in 
definierten Zeitpunkten (t) eines vorgegebenen Zeitbereichs und interpoliert 
die Abtastwerte durch eine stetige Funktion (f(t)), wobei die stetige Funktion 
(f(t)) durch eine Summe einer vorgegebenen Anzahl (n e N) von Wavelets 
(W) gebildet ist und wobei jedes Wavelet (W) dem Produkt eines Abtastwerts 

mit einer Spaltfunktion (^^) mit einer Gauss'schen Glockenkurve (e^ , a 

e R) entspricht. 



(Rg. 2) 




Fig. 1 (Stand der Technik) 
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